Spontaneous symmetry breaking in the linear sigma model 
at finite chemical potential: One-loop corrections 
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We investigate spontaneous symmetry breaking within the linear sigma model with the SU(2) x 
U(l) internal symmetry at finite chemical potential, which was suggested as a model for kaon 
condensation in the CFL phase of dense quark matter. One-loop corrections to the scalar field 
effective potential as well as its propagator are calculated. Particular attention is paid to the type- 
II Goldstone boson that appears in the Bose-Einstein condensed phase. Furthermore, we show that 
the type-I Goldstone boson — the superfluid phonon — is allowed to decay due to the nonlinearity 
of its dispersion relation at high momentum, and determine its decay width. 
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I. INTRODUCTION 

One of the most striking and general consequences of 
spontaneous symmetry breaking is the existence of gap- 
less modes, guaranteed by the celebrated Goldstone the- 
orem pi . The number and properties of the Goldstone 
bosons are crucial for the low-energy dynamics (and the 
low-tempcrature thermodynamics) of the system. In par- 
ticular, these soft modes play a major role in transport 
phenomena such as heat conductivity or viscosity. 

It has been known for a long time that the physics of 
spontaneous symmetry breaking in Lorentz-noninvariant 
systems may be quite intricate. The fundamental result 
in this respect was achieved by Nielsen and Chadha [!|: 
They showed in a nonperturbative manner that, under 
certain technical assumptions, the energy of the Gold- 
stone bosons stemming from the spontaneous symmetry 
breaking is, in the long-wavelength limit, proportional 
to some power of momentum. The Goldstone bosons are 
then classified as type-I, if this power is odd, and type-II, 
if it is even, respectively. The numbers of the Goldstone 
bosons of the respective types are related to the number 
of broken symmetry generators by the following inequal- 
ity, 

# type-I GBs + 2 x # type-II GBs 

> # broken generators. (1) 

Eq. shows that whenever there is a type-II Gold- 
stone boson, the number of Goldstone bosons may be 
smaller than the number of broken generators. A pro- 
found example of a system where this happens is pro- 
vided by the ferromagnet. In the past decade, however, 
several other systems with type-II Goldstone bosons have 
been studied, ranging from various high-density phases of 
QCD 0, H, H, and the relativistic nuclear ferroma gnet 
|3 to Bose-Einstein condensed atomic gases 0, flpt lTT| 
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and the general effect of relativistic vector condensation 

In our recent paper [l4j we investigated the Gold- 
stone boson counting in a particular class of Lorentz- 
noninvariant systems — the relativistic linear sigma 
model at finite chemical potential. Based on previous 
partial results 0, 0] we clarified the connection of Gold- 
stone boson counting with the possibility of nonzero den- 
sities of Noether charges and conjectured a general count- 
ing rule: Non-zero density of a commutator of two bro- 
ken generators implies one type-II Goldstone boson with 
a quadratic dispersion law. As a consequence, the in- 
equality QJ is saturated, up to exceptional cases like the 
phase transitions, where the phase velocity of a type-I 
Goldstone boson may vanish, thus making it an 'acci- 
dental' type-II one. Within the linear sigma model, we 
were able to give a detailed proof of these statements. 

Our analysis was, however, purely classical — we 
worked all the time at the tree level. The present pa- 
per is intended to fill this gap. There are several reasons 
that make the inclusion of the loop corrections necessary. 
First, it is only at the one-loop level that the \(4> (f>) 2 in- 
teraction plays a nontrivial role. Indeed, at the tree level, 
it has no effect on the spectrum in the normal phase, 
while in the Bose-Einstein condensed phase it merely 
serves to stabilize the scalar field potential. Second, the 
loop corrections may play a significant role, especially in 
the vicinity of the phase transition, where the classical 
analysis cannot be trusted. Third, it is well known that 
when the static part of the Lagrangian has higher symme- 
try than the full Lagrangian, the tree level is not sufficient 
to determine the spectrum (even qualitatively) correctly 
|16|. The quantum corrections are therefore necessary 
in order to check the saturation of the Nielsen-Chadha 
inequality Q. 

For simplicity, we do not follow the general symmetry- 
breaking pattern of Ref. . Instead, we analyze a par- 
ticular model — the linear sigma model with the global 
SU(2) x U(l) symmetry. This model has been used to 
describe kaon condensation in the Color-Flavor-Locked 
phase of dense QCD 0,0. It is also a special case of 
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a class of models of relativistic Bose-Einstein condensa- 
tion investigated by Andersen \vi\ . 

The plan of the paper is as follows. Using the man- 
ifestly invariant formalism of generating functionals, we 
determine in Sec. [H]the one-loop effective potential and 
numerically find its minimum as a function of the chem- 
ical potential. In Sec. II I II we calculate the one-loop cor- 
rection to the scalar field propagator. We demonstrate 
analytically the presence of gapless (Goldstone) poles and 
then compute numerically the corrections to the tree- 
level dispersion relations of the type-II Goldstone boson 
and its massive counterpart. In Sec. IIV| we calculate 
the decay rate of the superfluid phonon [18J. The paper 
is concluded with a summary and a discussion of the re- 
sults. Some necessary technical details are deferred to 
the appendices. 



II. EFFECTIVE POTENTIAL 



are two Goldstone bosons, one type-I and one type-II. 
Their low-energy dispersion relations read 



/ \x 2 — M 2 

s = Vv^ |p| 



and E 



respectively. In addition to the Goldstone bosons, 
there are two massive excitations with energy gaps 
v /2(3/x 2 - M 2 ) and 2/x. 



A. One-loop correction 

Let us now proceed to the one-loop calculation. In 
order to account for the breaking of the U(l) symme- 
try associated with the particle number (or strangeness, 
in the context of kaon condensation) we introduce the 
formal Nambu doublet 



We work with the model of Miransky and Shovkovy Q 
and Schaefer et al. 5] which is defined by the Lagrangian 



(2) 



[The index i refers to the doublet representation of the 
SU(2).] The one-loop effective action is then given by the 
textbook formula l2(M. 



Here is a complex doublet of the global SU(2) sym- 
metry and the chemical potential /i, associated with the 
global U(l) symmetry (particle number), is included in 
the covariant derivative, D^cf) — (<9 M — i8^fx)(j) [l^. This 
model has been analyzed, at the tree level, in detail in 
Refs. 0, so we only summarize the main results for 
later reference. 

When ji> M the static potential develops a nontrivial 
minimum and the scalar field condenses. As a result, 
the SU(2) x U(l) symmetry of the Lagrangian J2J breaks 
down to its U(l) subgroup [different from the original 
U(l)]. Associated with the three broken generators there 



r 1L [<Mt] ^[^R-logdetA- 1 (4) 

plus counterterms. Here S[<j>, <f>'] — J d A x C(x) is the clas- 
sical action and A is the tree-level matrix propagator 
defined by 



S 2 S 



Displaying explicitly the matrix structure in the Nambu 
space, it reads, in the momentum representation, 



J 



A^b) 



[(Po + m) 5 



-M 2 
-2A<A 



-2A0t0]£i 
Ut 



[{po - M) 2 



—2X(j>i(j)j 
M 2 - 2 



(5) 



r 



Here the classical field cf> is already assumed to be con- 
stant. There is no lack of generality in this requirement 
as long as the vacuum is translationally invariant. 

With the constant classical field <f> we may evaluate the 
one-loop effective potential as 



Vvl = V c \ - - 



logdet A _1 (fe) 



(6) 



2 J (2tt) 4 

plus counterterms, where V^i = (M 2 — fi 2 )cp* cf> + A((p* cp) 2 



In this form, the effective potential is manifestly 
SU(2) x U(l) invariant. For detailed calculations it is, 
however, more convenient to fix the direction of <f> in the 
SU(2) 'flavor' space. As usual, we set 



-( ° 



(7) 



With this choice, A 1 becomes diagonal in the flavor 
space and may easily be inverted, yielding 
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An(p) 
A 22 (p) 



1 ( (pa - ii) 2 - p 2 - M 2 - 2\v 2 Xv 2 

(Po + M) 2 P 2 



bg-^(p)][^-^ 2 (p) 



Aw 2 
I — 



M 2 - 2Av 2 



(8) 



where 

E 2 ± {p) = p 2 + n 2 + M 2 + 2Xv 2 

± v/4/i 2 (p 2 — n 2 + M 2 + Xv 2 ) + [2y? + Xv 2 ) 2 . 

By minimizing the classical potential V c \ we find that, at 
tree level, Aw 2 = /.i 2 — M 2 , which yields the dispersion 
relations 



E 2 = p 2 + 3/i 2 - M 2 ± P 2 + (3^ 2 - ^ 2 ) 2 , (9) 

in accord with the result of Refs. (HIH- 



B. Renormalization 

The loop integral in Eq. © is divergent. To renormal- 
ize it, we add to the effective potential the counterterm 



counterterm 



Sv + SM 2 ^<p + 5A((/> t 



This term generates, by means of its functional deriva- 
tives, also the counterterm for the scalar propagator. The 
constant Sv is irrelevant for our purposes as it just fixes 
the offset of the effective potential at </> = 0. The coun- 
terterms in the effective potential are sufficient to absorb 
all infinities for, as is well known, there is no wave func- 
tion renormalization in the scalar 4 theory at one loop 

Actually, we shall not need to evaluate the effective 
potential at all. What we are interested in is just the 
vacuum expectation value of the scalar field, or v. We 
shall be therefore solving the 'gap equation' dVy^jdv = 
0. Now since the effective potential is a function of v 2 , 
there is always the solution v = 0. In the Bose-Einstein 
condensed phase this is, however, not the only solution 
and not the global minimum. The nontrivial solution is 
found by means of dVy^jdv 2 = 0. Formulas JBJ and JHJ) 
give 



dV, 



1L 



1 



1 



dv 2 



= -(M + SM 2 - p ) + -(A + SX) 



2AA?f (k) 



g v 2 det A 22 



l(k) 



det A^(k) 



(10) 



The momentum integral in Eq. I] 10(1 is only formal and 
implicitly involves a convenient regularization. In deriv- 
ing Eq. (|10|) we already assumed that the regularization 



allows to change the sign of the integration variable so 
that the contribution of the particles circulating in the 
loop is the same as that of antiparticles. For the later 
calculation of the scalar field propagator it will also be 
useful to be able to shift the integration variable. In con- 
crete computations we shall use the analytic integration 
over energy (frequency) in combination with dimensional 
regularization around d — 3 space dimensions |17| . 

Eq. HlOfl may be rewritten in a particularly useful way. 
With a decent use of Eq. we find (for the sake of 
brevity, we temporarily drop out the arguments of the 
propagators as well as the integration measure) 



^detA- = ^(A-« 



A- 



a; 



"*A 



= -2A ( A 



A: 



-a; 



so that 



dv- 



■ det A 



22 



det A^ 1 



= 2A 



Atl* + Atf + A 



2A / 2A 



22 ' '-'22 

t 



■*22 



Eq. (|10fl thus acquires the form 
^ 2 - (A/ 2 + SM 2 ) - (A + SX)v 2 



2iX 



d*k 



Atf (k) + 2Al* J (k) + Ali(k) (11) 



Even though this result has been derived by a not very 
transparent manipulation, it could have been expected: 
It is equivalent, up to an overall factor, to the require- 
ment that the one-particle-irreducible tadpole contribu- 
tions to the vacuum expectation value of <f>2 vanish. For 
the diagrammatic representation of the individual terms 
see Appendix iBl 

For the concrete calculations, we adopted the following 
renormalization conditions. The scalar mass was renor- 
malized by subtracting the whole one-loop correction at 
[i = and v = [see also Eq. (|T3|) later], i.e., 



SM 2 



-3AA 



3-d 



d d k 



1 



{2n) d Vk 2 + M 2 ' 



(12) 



where the integral is regularized by minimal subtraction 
in d dimensions and A is the renormalization scale. This 
procedure guarantees that the parameter M keeps its 



4 



interpretation as the physical (pole) mass of the scalar 
at zero chemical potential. 

The coupling constant A was renormalized by modified 
minimal subtraction (MS) so that 



0.01 



sx = 



4tt 2 



- — 7 + log 47T 



where e = 3 — d is the expansion parameter of dimensional 
rcgularization 28]. (The pole does not depend on the 
chemical potential and hence could be extracted from 
the analytically calculable loop integral at [i = 0.) The 
renormalization scale A was set equal to the characteristic 
scale of the system: we chose A = M in the normal phase 
and A = /i in the Bosc-Einstein condensed phase. Since 
the phase transition occurs at fi = M (see the following 
discussion) , the subtraction point is a connected function 
of the chemical potential. 

For sake of numerical solution the gap equation 
was rewritten, with the help of Eq. (JSJ and upon the 
integration over frequencies, as 



v 2 = 



fi 2 - M 2 



1 



d 3 k 

(2tt) 3 I Vk 2 + AP + Xv 2 

k 2 + M 2 + |Ai; 2 -/! 2 i 



E + (k)+E„(k) 



1 



£+(k)£_(k) 

SM 2 + SX i 



A 



(13) 



This equation was solved iteratively with the initial 

ansatz given by the tree-level value, vq = yj M 2 ~ M , To 

implement the MS scheme within the numerical compu- 
tation, we subtracted the pole part in the form 



« 3A2 l A2 



3X Z A 



2 a 3-d 



d d k 



1 



(27T) d (k 2 + M 2 ) 3 /2 



The integral here reproduces the pole and is, in fact, 
equal to the full one-loop correction to A at /x = and v — 
0. Upon the subtraction of SM 2 and SX, the integral in 
(|13|l is rendered finite and was evaluated approximately 
with a simple cutoff. 

The numerical results are displayed in Fig. ^ We em- 
phasize the fact that the phase transition to the Bose- 
Einstein condensed phase happens at the chemical poten- 
tial equal to the renormalized scalar mass. This is of 
course not surprising from the physical point of view and 
in fact is an important check of consistency of our calcu- 
lations. 

This conclusion is also easily proved analytically. Just 
set v = and then immediately observe that the right 
hand side of Eq. I|ll|) is exactly canceled by the coun- 
terterm SM 2 i.e., v = is a solution of Eq. (jll|) exactly 
for n = M. 

Finally, note that the convergence of our iterative so- 
lution of the gap equation guarantees that the found so- 
lution is, as required, a (at least local) minimum of the 
effective potential. 




FIG. 1: Relative increase of the vacuum expectation value 

v with respect to the tree-level value vo = \J ^ ~ X RI - ■ The 
numerical data were obtained with the renormalized coupling 
A set to 1. 



III. SCALAR PROPAGATOR 

The one-loop-corrected propagator V is obtained as 
a second functional derivative of the one-loop effective 
action 10} , 



S 2 T 



1L 



S^S^iy) 



After differentiating the logarithm of the determinant we 
arrive at 



■Tr 



A 



S 2 A- 



. SA- 1 . 8 A- 1 
- A — 7 — A- 



plus counterterms. Here 'Tr' denotes the full functional 
trace, over both spacetime and internal degrees of free- 
dom. 

The derivatives of the inverse bare propagator A -1 are 
nothing but the cubic and the quartic interaction vertex. 
They are local and may be represented by the coupling 
matrices T« and defined 



tot Or) 
* 2 A fc 7(M) 



= T$S i (s-t)S i (s-x) 



Q t ti ) 5\s-t)5\s-x)8\t-y). 



The explicit formulas for the matrices and are 
listed in the Appendix lAl With this notation, the inverse 
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propagator acquires the form 



13 

2 / To^U 



j£k_ 

(2tt)' 



tr 



A(k)Q {ij) - A(k)T^A(p + k)T^ 



(14) 



plus counterterms. The symbol 'tr' here refers to a trace 
over internal degrees of freedom only, that is, over the 
flavor and Nambu spaces. 



Note that the loop correction in Eq. I|15fl i s 
momentum-independent, as usual in a </> 4 theory. The 
renormalization is therefore trivial and corresponds to 
a mere redefinition of the mass. 

Since in dimensional regularization we may freely shift 
the integration variable, the loop integral in Eq. 11151) 
is //-independent and thus is completely canceled by the 
mass counterterm (|12[l . Therefore, in the normal phase, 
the scalar field propagator gets no correction and the 
mass spectrum simply consists of two doubly degenerate 
levels atM±/i 



iimpr 

m 



A. Normal phase 

When the chemical potential is small enough, the sys- 
tem is in the normal phase: v — and the symmetry 
is not spontaneously broken. In this phase the one-loop 
correction to the propagator is particularly simple. The 
bare propagator A is diagonal in both the flavor and the 
Nambu space, the only nonzero components being 



M 



2l-l 



A^(p) = [( P0 - M ) 2 -p 2 -M 2 ]- 1 . 

Since in the normal phase there is no cubic coupling 
TW, the final formula for the one-loop propagator is eas- 
ily found to be 



(P) = [(Po + V? - P 2 - (M 2 + 6M 2 )} Si, 
d A k 1 



6i\5i 



(2tt) 4 (fco + M) 2 - k 2 - M 2 ■ 



(15) 



The other nonzero component of the propagator, D^t^, 
is related to Eq. (|15|) by the first of the identities, 



V, 



(p) = = vU-p)- 



(16) 



which follow from the very definition of the propagator. 



B. Bose—Einstein condensed phase 

When the scalar condenses i.e., v > 0, we have to 
work with the full formulas ©■ The matrix propagator 
T> has altogether 16 components (four flavor times four 
Nambu). Fortunately, several of them are actually zero 
due to the clever choice of the vacuum, Eq. (0. With this 
choice, the unbroken flavor U(1)q symmetry is generated 
by the matrix Q = i(l + t 3 ). This means that the upper 
flavor component of <p, <j>\, carries the unbroken charge 
while the lower component </>2 does not. The charge con- 
servation then immediately implies that the propagator 
is diagonal in the flavor space, i.e., T>±2(p) — ^21 (p) = 0. 
Moreover, the propagator of <fii is diagonal in the Nambu 

space, that is, T>ff(p) = Vf^ (p) = 0. [All these rela- 
tions could, of course, also be demonstrated explicitly by 
a proper analysis of the formula (|14|) .] 

The remaining nonzero components of the propagator 
are strongly constrained by Eq. (|16J) so that only three 
of them are independent: T>ff , Djj and Uja ■ 



1. Propagator of 4>i 

The only independent component of T>u is T>ff 
A straightforward application of Eq. (|Tlf> yields 



V 



11 
2i\ 



(p) = [p a + (if 
(2^ 



(M 2 + SM 2 ) 



2Aff (fc) + Afl (k) 



2i\ 2 v 2 



(A + 5X)v 2 
d A k 



Aff (p + k) AH (k) + Aff (A;) + A^ (k) + A£»(fc) 



(17) 



The first integral represents the tadpoles and is 
momentum-independent. A qualitative difference in 
comparison with the normal phase comes in the second 
integral, which contains the contribution of the cubic ver- 
tices. The diagrammatic representation of all terms in 



the propagator is given in Appendix iBl 

The propagator vff (p) is expected to have a mass- 
less pole corresponding to the Goldstone boson. Indeed, 
a short glance at Eq. JHJ reveals that at tree level, when 

Xv 2 = /i 2 — M 2 , there is a gapless pole in Aff (p). The 
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field 4>i annihilates a particle with dispersion relation 
E = \/p 2 + /i 2 — fx, i.e., a type-II Goldstone boson. The 
other excitation, annihilated by (j>\ , with a tree-level gap 

2/i is manifested as a pole in T>f 1 ^(p). 

The analytic proof of the existence of a massless pole 
in T>ff (p) is provided in Appendix IC II Here we just 
observe that once the relation T>^^ (0) = is proved, 

we may subtract it from T> (p) to obtain a convenient 
expression 



V 



n ' (P) = Po + 2 t i Po 



2i\ 2 v 2 



d 4 k 



Aff (p + fc)-Aff(fc) 



x A** (A) + A#(fc) + A** (fc) 



A^(fc) 



(18) 



Note that the propagator no longer depends explicitly on 
the renormalization counterterms. Also, the momentum- 
independent tadpole graphs disappeared and the remain- 
ing integral is finite and perfectly well defined even when 
the regularization is removed. 

In order to establish the radiative corrections to the 
type-II Goldstone boson dispersion relation, we per- 
formed the integration over frequency and expanded the 
result in powers of the external momentum p, up to or- 
der po and p 2 , respectively. The final formula is rather 
cumbersome, so we do not write it out and instead re- 
port just the result of the numerical integration of the 



remaining momentum integral. 

At the leading order of the momentum power expan- 
sion, the inverse propagator reads 



V 



li 



(p) = 2[ip (l + Z 1 )-p 2 (l + Z 2 ). 



The coefficients Z\ , Z 2 determine the renormalization of 
both the wave function and the dispersion relation of 
the Goldstone boson. Their dependence on the chemical 
potential is displayed in Fig. [21 

Now we turn our attention to the massive partner of 
the type-II Goldstone boson. Analogously to Ref. 
we shall calculate just the leading order of its disper- 
sion relation, i.e., its gap. The corresponding propa- 
gator, , is obtained from Eq. I|18fl y i a the rela- 
tions (|16fl . Next we set p = 0. However, the explicit 
form of the tree- level propagators, Eq. (JSJ, reveals that 

A?f (fco — 2//, k) = Aji (fco,k) so that the difference of 
propagators in the loop correction to T>ff (—po, 0) ex- 
actly vanishes at po = 2/i, see Eq. (|18H . As a result, the 
renormalized gap is equal to its tree-level value 2\x and 
receives no radiative corrections. 



2. Propagator of <j>2 

In this case, Eq. I|14(l yields, after a reasonable appli- 
cation of the rules lfRI)l . 



v^t"f\) = (p + »f 



(M 2 + SM 2 ) - 2(A + 6X)v 2 - 2i\ 



d 4 k 



Aff(k) + 2Alf(k) 



2iX 2 v 2 



(27T) 



4 A?f {k)A\f (p + k) 



8iAV 



sin* 

(2tt)4 1 22 



(p + k) 



Atf (k) + A||(fe) + Af/ (k) + ~A^(fc) 



+ Ait( P + k)Ai 2 Uk)\, (19) 



V^( P ) = -(A + SX)v 2 2 l X I ^jA^fc) - 2iAV 



d 4 k 



8iX 2 v 2 



d 4 k 



(27T) 



Att(p + k) 



' rk Atf{k) Aff {p + k) 



(27T) 



4t\k) + -Att{k) + Af/{k) 



\k) 



+ Att (p + k)Atl (k) . (20) 



r 



Here the propagator is not diagonal in the Nambu 
space and the situation is thus more complicated. We 



Fortunately, we need not calculate the determinant ex- 
plicitly, at least to prove the existence of a Goldstone 



cannot look for the zeros of the individual components of boson. Recalling the relations (|16fl , the determinant at 



the inverse propagator. Instead, we have to solve the con- 
dition detT>22 (p) — 0, which seems a rather formidable 
task in view of Eqs. JT5J and (f2U|) . 



zero momentum simplifies to 



detD 2 - 2 1 (0)= fe^O)! - V-^iO) 
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FIG. 2: Renormalization constants of the type-II Goldstone 
boson propagator. The numerical results indicate that their 
dependence on the chemical potential is linear near the phase 
transition and saturates at large fi. The data were obtained 
with A = 1. 



The secular equation thus reduces to 



In view of Eq. (0 we see that the existence of a Goldstone 
boson is proved once we verify the relation 



(21) 



A detailed proof of this statement is provided in Ap- 
pendix E21 

The analysis of Eqs. (fTT)|l and (|2T)|l would yield the ra- 
diative correction to the type-I Goldstone boson disper- 
sion relation, i.e., the speed of sound in the Bose-Einstein 
condensate. We do not perform the detailed calculation 
here since our main concern was the type-II Goldstone 
boson studied above. Moreover, such a straightforward 
approach is unnecessarily complicated because the one- 
loop correction to the speed of sound has been deter- 
mined by Andersen using a different method |17| . In the 
following, we shall instead concentrate on another spe- 
cific property of the superfluid phonon in our system — 
its decay. 



IV. PHONON DECAY 

Upon inspection of Eq. ©, the full tree- level dis- 
persion relation of the type-I Goldstone boson (or, the 
phonon) turns out to be convex. This fact has far- 
reaching consequences for it kinematically allows the 
phonon to decay. (A similar, convex, dispersion law is 
found e.g. for the superfluid phonon in the Bose-Einstein 
condensed weakly interacting Bose gas [2lJ.) This in turn 



strongly affects the transport properties of the system 
such as the shear viscosity |2^|. 

Here we calculate the two-particle decay rate of the 
phonon at the leading order of perturbation theory. Such 
a quantity could in principle be extracted, via the opti- 
cal theorem, from the imaginary part of the one-loop 
propagator calculated in Sec. IIII B 21 This would, how- 
ever, be rather complicated for two reasons. First, the 
phonon is, in the propagator of 02, entangled with the 
massive 'radial' mode. Second, we would have to deal 
with the full momentum dependence of the propagator; 
the low-momentum limit studied before would not be suf- 
ficient for our purpose. Since the upcoming calculation 
of the phonon decay rate involves a few rather nonstan- 
dard pieces, we regard it instructive to perform it step 
by step, using ordinary perturbation theory. 

Before starting the calculation we briefly comment why 
we have chosen to study such a particular process. First, 
as has already been mentioned above, the low-energy 
properties of the Goldstone bosons are essential for the 
thermodynamics of the system at low temperature. Two 
factors are determinant for the possibility of the decay of 
the Goldstone bosons: The conservation of the unbroken 
U(1)q charge and phase space (kinematic) restrictions. 
The type-II Goldstone boson, G, carries the unbroken 
charge. Since at low energy it can only decay to gapless 
particles, the only processes allowed by charge conser- 
vation arc G — > Gtt, G — ► G7T7T, etc. (Here and in the 
following, 7T denotes the type-I Goldstone boson i.e., the 
phonon.) Taking into account the dispersion relations 
©, these decay modes, however, fail to satisfy the en- 
ergy and momentum conservation simultaneously. Thus, 
the type-II Goldstone boson is stable at low momentum. 

On the other hand, the phonon may decay, by the pro- 
cesses 7r — > 7T7T, 7r — > tttttt etc. Only the two processes ex- 
plicitly written here receive nonzero contribution at the 
first order of perturbation theory. As will become clear 
later, the three-particle decay is suppressed, at least in 
certain ranges of momentum and chemical potential. We 
shall therefore calculate just the decay rate for the pro- 
cess 7r — > 7T7T and comment on the relevance of our result 
afterwards. 



A. Vacuum transition amplitudes 

The decay rate cannot be straightforwardly deduced 
from the Lagrangian @ because it is expressed in terms 
of the scalar field (j> rather than the physical modes n, G 
and their massive counterparts. The problem is that the 
kinetic term in the Lagrangian cannot be diagonalized at 
the level of fields. Indeed, Eq. JHJ shows that the di- 
agonalizing unitary transformation is necessarily energy- 
dependent. We therefore have to calculate the three- 
point Green's function of <f) and use the LSZ reduction 
formula to extract the physical scattering amplitude. 

Alternatively, we may directly calculate the S-matrix; 
the couplings of cj> to the physical states follow from 
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the Kallen-Lehmann representation of the propagator. 
Let A(x) and B(x) be two local bosonic operators. 
The general form of the Kallen-Lehmann representa- 
tion of the time-ordered Green's function, T> AB (x,y) — 
—i(0\T{A(x)B(y)}\0), reads, in the momentum represen- 
tation, 



(0|A(0)|n,p)<n,p|fl(0)|0) 



lu — E(p) + ie 

fl(0)|n,-p)(n,-p|A(0)|0) 
lu + E(p) — ie 



(22) 



where E(p) is the energy of the Hamiltonian eigenstate 
|n,p). The summation index n is discrete for one- 
particle intermediate states, and continuous for multi- 
particle states. Also, the states are assumed to be nor- 
malized as (n, p\m, k) = S mn 6 3 (p — k). Note that we use 
the general Kallen-Lehmann representation 122(1 rather 
than the covariant form usual in literature |2(l |. 



The poles in the propagator of 4> 2 correspond to the 
dispersion relations given in Eq. The contribution 

of the massless pole at lu = E- (p) — ie is 



A| 2 ole (p, 



2E_ (p) Ei (p) - Ei (p) uj - E- (p) + ie 



[E-(p) - /i] 2 - p 2 - M 2 - 2Xv 2 
Xv 2 



Xv 2 

[£L(p) + a*] 2 - P 2 - M 2 - 2Xv 2 



This expression immediately yields the transition amplitudes 

1 1 1 



<O|0 2 (O)|tt(p)) = 



= {-[£L(p) - fi} 2 +p 2 + M 2 + 2Aw 2 } 1/2 , 



(2tt)3/2 j2EUv) ^E 2 + {p)-EUp) 

(0|4(0)K(p)) = - 1 1 ={-[E-(p) + fA 2 + P 2 + M 2 + 2Xv 2 } 

(2tt)3/ 2 y/2E..(p) ^E 2 + (p)-E 2 _(p) 1 J 

I 



1/2 



(23) 



The relative sign of the amplitudes follows from the non- 
diagonal part of the propagator. Note that the ampli- 
tudes are particularly simple in the limit v — > (i.e., at 



the phase boundary) where E±(p) = \[p 2 
e p ± [i. It follows that, in this limit, 



M 2 ±n 



(O|0 2 (O)|tt( P )) 



(2tt)3/2^ 



(0|4(0)|7r(p)) = 0. 



Thus the result typical for Lorentz-invariant theories is 
recovered in the normal phase. 

At the first order of perturbation theory, the ampli- 
tude for the two-particle decay of ir follows from the La- 
grangian upon shifting 2 by its vacuum expectation 
value, v/y/2. The relevant cubic terms read 

-V2Xv(<f>l<f> 2 4>2 + 4>thh)- 

These two terms generate three diagrams each, contribut- 
ing to the given decay process, that differ just by a per- 
mutation of the lines. The result for the S- matrix ele- 
ment is 



S = -{2TT) A 8 A ( P -k-q)iM 



I 



[(2,r) 



3/21 



where the scattering amplitude M. is given by 



-iM- 



= -2V2iXv 



[(2703/* 

'<7r(k)|4|0)(7r(q)|4|0)(0|^|7r(p)) 

(7r(k)|4|0)(7r(q)|^|0)(0|4|7r(p)) 
(7r(k)|0 2 |O)(7r(q)|4|O)(O|4|7r(p)) 
(7r(k)|0 2 |O)(7r(q)|0 2 |O)(O|4|7r(p)) 
(7r(k)|0 2 |O)(7r(q)|4|O)(O|0 2 |7r(p)) 
+<7r(k)|4|0)<7r(q)|^|0)<0|^|7r(p)) 



(25) 



For legibility we omitted the arguments '(0)' of the field 
operators. Also, the initial state is labeled by the mo- 
mentum p, the momenta of the products are k, q. 



B. Kinematics 

Our analysis of the two-particle decay is further com- 
plicated by the fact that the Lorentz invariance is broken 
by the chemical potential. We therefore have to evalu- 
ate the decay rate as a function of the momentum of the 
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decaying particle. After the common trick of squaring 
the S-matrix in finite volume and performing the trivial 
three-momentum integration we arrive at the formula for 
the decay rate, 



r = 



1 1 

24^ 



J k 2 d 3 k\M\ 2 5(E-(p) - E_(k) - E_{q)). 



The factor ^ arises from the fact that the two par- 
ticles in the final state are identical. The ^-function 
fixes the moduli of the momenta k and q as a func- 
tion of the scattering angle 9, which measures the de- 
flection of k from the direction of p. The |q| is given by 
|q| = ^/p 2 + k 2 — 2|p||k| cos# and |k| is then fixed by 
energy conservation. Due to the form of the dispersion, 
Eq. © , the solution of the energy conservation condition 
is, however, rather involved, so we work it out analyti- 
cally only in a special case (see Sec. IIV C|l . In the general 
case, the decay rate is computed numerically. 

Using the chain rule to differentiate -E-(q) in the 6- 
function, we arrive at the final formula for the decay rate, 



1 

47T 



sin6>d(9k 2 (<9) 



E'_{k(8)) + E'_{q(8)) 



|k(0)|-|p|cos0 



(26) 



where m is the maximum angle to which n can de- 
cay at a given energy. Fig. |21 displays the expression 
— S_(p) + -E-(k) + -E-(q) as a function of |k|, for sev- 
eral values of the angle 9. The critical angle 8 m is ap- 
parently determined by the condition that the deriva- 
tive of the argument of the energy-conserving (5-function 
[i.e., the expression in the denominator in Eq. (|26ll ] 
vanishes at k = 0. Since for k = both energy and 
momentum conservation require q = p, this reduces to 
cos# m = E'_{Q) / E'_(p) . Using explicitly the dispersion 
relation ©, we have 



cos 9„ 



— M 2 



3/i 2 - M 2 



2n 2 



Wp 2 + (3/i 2 - M 2 ) 2 



(27) 



Formula (|27|l admits several particularly simple lim- 
its. For |p| — * we find 9 m = 0. This is easily un- 
derstandable: At low momentum the phonon dispersion 
relation is linear and the decay thus kinematically for- 
bidden. For |p| — * oo, on the other hand, we have 

3 ^ 2 ^-M 2 • Finally, for /x — > M+ at fixed 
momentum i.e., at the phase transition, 9 m = ir/2. This 
shows that one has to be careful when performing these 
limits. Near the phase transition there may be two small 
scales, |p| and /i — M, and the result depends on their 
hierarchy. 
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FIG. 3: Sample plot of the expression — -E-(p) + E-(k) + 
_E_(q) (in units of M) as a function of |k|, for [i = 2M and 
p = M. The curves are labeled by the value of the angle 9 
in degrees. The critical angle for the chosen set of parameters 
is 8 m « 26.3 deg. 



C. Decay rate 

The phonon decay rate may be determined analytically 
in the limit v — > 0. Strictly speaking, the following cal- 
culation is justified only for p 2 ^ Xv 2 . However, later 
the analytic result will be checked by a direct numerical 
integration of the full formula l|26|) . For the time being, 
we simply set v — everywhere except, of course, the 
overall factor v in the amplitude M.. 

Eq. (|24f> shows that only the first term survives in Eq. 
PSfl so that 



M = 2V2XV- 



1 



1 



1 



Setting \i = M in Eq. © gives ^(p) — e p ± M, hence 
E'_(p) = |p|/e p and Eq. lj2l)|) reduces to 



A 2 2 «r/2 

T = / sin9d9 

47re P Jo 

X 



k 2 (9) 



H°)\( e k(e) +e q (0)) - |p|ek(e)Cos6>' 



(28) 



With the simplified dispersion relation, the energy and 
momentum conservation may easily be solved to yield 



2M|p|£ + (p)co^ E _ {k) 
El(p)-p 2 cos 2 9' v ' 



2Mp 2 cos 2 9 
E\{p) — p 2 cos 2 



The angular integration in Eq. I|28|) is then elementary 
and we arrive at the final formula for the decay rate, 



2„,2 



F 



\ 2 v 



4tt e p (e p + M)' 



(29) 
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FIG. 4: Decay rate for the process n — > nn as a function of 
momentum (solid lines). The numerical data were obtained 
with A = 1. The curves are labeled by the dimensionless pa- 
rameter (x/M — 1. The dashed line corresponds to the analytic 
result, Eq. I|290 , The decay rates for different values of the 
chemical potential are scaled by the factor Xv 2 so that their 
convergence to the analytic result is comparable. 

Such a dependence of the decay rate on the momentum 
of the initial phonon might have been expected: The de- 
cay rate tends to zero at both small and large momenta, 
where the dispersion relation becomes linear. A compar- 
ison with the exact numerical integration of Eq. (j26(l is 
displayed in Fig. 0] 

V. CONCLUSIONS 

In this paper we studied spontaneous symmetry break- 
ing by Bose-Einstein condensation in scalar tfi 4 theory 
with two charged scalar flavors. Our primary goal was 
to refine our previous results |14| by including radiative 
corrections. 

Within the one-loop calculation we verified that the 
Bose-Einstein condensation sets on at the chemical po- 
tential equal to the renormalizcd scalar mass (at ji = 0). 
This assertion, though perfectly physical, may seem sur- 
prising within the approximation used since, as is well 
known 20] , in massless <fi theory the one-loop correction 
draws the scalar field expectation value off the origin to 
a nonzero value. The spontaneous symmetry breaking 
thus achieved, however, eventually turns out to be an ar- 
tifact of the one-loop effective potential and disappears 
upon the renormalization group improvement. 

Fortunately, such a fake vacuum apparently does not 
occur in our case for the phase transition point is exactly 
where it should be. This gives us hope that the renor- 
malization group improvement, however important for 
critical phenomena, should not change the conclusions of 
this paper qualitatively. 



In Sec. IIII B II we determined the radiative corrections 
to the dispersion relations of the type-II Goldstone bo- 
son and its massive counterpart. The numerical results 
suggest that the dependence of the renormalization con- 
stants Z\ and Zi on the chemical potential saturates at 
large fi. Since in this region the two massive modes in 
the model become very heavy, the asymptotic behavior 
of the type-II Goldstone boson dispersion relation should 
be governed by a low-energy effective theory for the Gold- 
stone bosons, presumably analogous to that for the non- 
relativistic ferromagnet 

Furthermore, we showed that the massive counterpart 
of the type-II Goldstone boson has gap 2/x, which re- 
ceives no radiative corrections. This assertion has a sim- 
ple physical interpretation. Recall that /i is, originally, 
a chemical potential associated with the U(l) symmetry 
of the Lagrangian Once this symmetry is sponta- 
neously broken, /i may be reinterpreted as a chemical 
potential of the unbroken U(1)q. Since 4>\ and <fi\ carry 
the charges Q = ±1, the corresponding excitations differ 
by 2[i in energy. 

In a sense, this is a remnant of the Goldstone theo- 
rem: It guarantees nonrenormalization of the Goldstone 
boson mass to all orders in the loop expansion. Once the 
chemical potential is switched on, the excitation annihi- 
lated by can no longer be massless, but the Goldstone 
theorem still assures that its gap does not renormalize. 

In Sec. IIVI we determined the two-particle decay rate 
of the superfluid phonon. Such a quantity as well as the 
very possibility of the phonon decay may be important 
for the transport phenomena in the Color-Flavor-Lockcd 
quark matter with a kaon condensate p^ . As shown in 
Refs. 0i EMEU j the CFL phase with a meson conden- 
sate appears to be energetically preferred to the simple 
CFL phase, thus making our results of a possible rele- 
vance for the phenomenology of compact stars. 

To summarize, in this paper we investigated sponta- 
neous symmetry breaking within the linear sigma model 
at the one-loop level. We verified the saturation of the 
Nielsen-Chadha inequality and hence also the count- 
ing rule proposed (and proved at the tree level) in our pre- 
vious work UJ. We thus fulfilled one part of the program 
outlined in Ref. |14|. The other, more ambitious part is 
to extend the results achieved so far to other relativistic 
systems with finite chemical potential. Our recent work 
shows that at least the type-II Goldstone bosons may be 
described, to a large extent, in a model-independent way. 
These results will be reported elsewhere. 



APPENDIX A: CUBIC AND QUARTIC 
COUPLING MATRICES 

In this appendix the explicit form of the coupling ma- 
trices T« and is given. The indices <fi and (f>^ refer 
to the upper and lower components of the Nambu dou- 
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blet, respectively. 



rp(i) _ 
± kl4 > — 

rp(i) 

± kl<f>i 

m 

m _ 



Q 
Q 



Q (lj) 



-2\(5u(j)k + 8u4>i) 
■ 5m4) 



-2X(5 U (t>l - 

-2A(<Wj 




-2X(Sj k 5 l i + SkiSij) 






—2X(SikSji + SjkSu 

-2X(6 ik Sji + 5jk5u) 


-2X(Sik5ji + SkAj) 






—2X{8ik4>i + 5ki4>i) 
-2X(5 a (j> k + 8 ik ct>i) 
-2X(5 U (l>l + S k i4) 


-2X(SikSji + SuSij) 








-2X{5 jk 5ii + SkiSij) 



We start with the interaction vertices. These all come 
from the — X((jr (j)) 2 term in the Lagrangian. Upon shift- 
ing 02 by v/^/2, the interaction yields the following ver- 
tices, 



V 
4 



y 

\ 
/ 

I 

\ 

4 



= -4iA, 



/ 



-2V2iXv, 



= -2iX, 

-V2iXv, 
'V / 

a 'v 



AiX. 



APPENDIX B: DIAGRAMMATIC 
REPRESENTATION 

In the following, we summarize the diagrammatic rep- 
resentation of the individual terms in the effective poten- 
tial as well as the propagators. The flavor components 
<pi and fa are denoted by oriented dashed and dotted 
lines, respectively. Since the D22 propagator is not di- 
agonal in the Nambu space, it has to carry two arrows, 
specifying the operators it connects at its ends. The bare 
propagators (JSJ) are represented by the following lines. 



A<W> T . 
^11 

^22 

^22 



Ml 



One-particle irreducible contributions to the vacuum 
expectation value of the scalar field 2 > see Eq. (jllf) : 



We have explicitly indicated the symmetry factor \ at the 
last graph. Note that in our formalism all such factors are 
automatically generated by the functional differentiation 
of the effective action. 

Next is the diagrammatic representation of the propa- 
gators. The loop graphs are in one-to-one correspondence 
with the terms in the formulas {T7|) , and JSUJl . 



u \\ — ^11 



4' \ 



4 > 



v- 



= a: 



1 4 



4' \ 



V / 



4 \ 



t > 



u 21 



\ f 



4 \ 
V 4 



V 4 



4' > 

v V 



V 4 



4 \ 

V. > 
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APPENDIX C: PROOF OF THE EXISTENCE OF GAPLESS POLES 



Here we provide a detail proof of the existence of massless poles in the propagators T>\\ and X>22- The key ingredients 
are the formulas l|17|) , l|19fl and (|20[1 , as well as the gap equation l|ll|l . 



1. Pole in ©11 



We are going to prove that Z? n (0) = 0. Eq, 117|) together with the expressions © for the bare propagators give 



£>22 100t (O) _ (j, 2 -(M 2 +5M 2 )-(X+8X)v 2 -4iX 



d 4 k 



(2tt) 4 (fc + m) 2 - 4 - A« 2 



-2iX 



d A k (fco — /i) 2 



2Xv 2 



(2tt)< [ki-Eumki-Eun 



2iX 2 v 2 



d A k 



Xv 2 } + {(k - nf 



Xv 2 



(2tt)4 (fc + m) 2 -e 2 k - Xv 2 



[k 2 -E 2 + (k)}[k 2 o~E 2 (k)} 



(CI) 



where we abbreviate ej. = Vk 2 + M 2 . Next we expand the gap equation Ijlll) , using the explicit form of the propagators 



■i 2 - (M 2 + 5M 2 ) - (A + 5X)v 2 = 2iX [ 

J (2^) 4 (fco + /i) 



,V 11 A J {2^(k + l i) 2 -ei-Xv 2 

f d*k (k + fi) 2 — 4 — 2Xv 2 f d*k 1 

J (2tt)4 [k 2 E 2 + (k)][k 2 E 2 (k)] + J (2tt)4 [fcg - E\ (k)][fc 2 - E 2 _ (k)] 

stitute for ^ 2 - (Af 2 + SM 2 ) - (X + SX)v 2 in Eq. l(Clj). which finally yields 

2Xv 2 |\ Aw 2 



" \ / LU "I" \ / J I. U \/J 

This expression is used to substitute for /i 2 — (M 2 + SM 2 ) — (A + <5A)v 
i^t,„, f d 4 k ( I (fcp - /i) 2 - c 



11 V ; J (2tt) 4 \ (fc + m) 

Using the identity 



Xv 2 [fc 2 -£ 2 (k)][fc 2 -£ 2 (k)] (fc + /i) 2 - e 2 - 2Xv 



det A-^fc) = [fc 2 - £ 2 (k)][fc 2 - E 2 _ (k)] = [(fc + /i) 2 - 4 - 2At; 2 ][(fc - m) 2 - e£ - 2A V 2 ] - (At; 2 ) 2 , 

'and in the last expression vanishes. 



A 22 - j '■ ' ' * 

which follows from Eq. (jSJ) , it is now straightforward to show that the integr 



2. Pole in E> 22 



In this case we are going to prove Eq. (|21l) . We shall use a more symbolic method, without referring to the explicit 
form of the bare propagators |JS}. Upon a comparison of Eqs. H19f) and 12U|) we can see that several terms immediately 
cancel out. After a mild rearrangement we are to prove the identity 

f _ (M 2 + (5M 2 ) _ (A + sx y = 2lX J + 2A ^t _ A ^ + 4a 2 w 2 J ( A *f a£* _ A#A#) . (C2) 

dct A22 

For the sake of legibility we dropped out the integration measure and arguments of the propagators. Also, we used of 
the fact that = A^. 

On the right hand side of Eq. (|C2() . we rewrite Aw 2 detA22 as Aw 2 / det A^ 1 = A%2 and immediately arrive at 
an expression identical to the right hand side of Eq. . Eq. (|C2(I is thus proved and hence also the existence of 
a massless pole in the propagator 2?22- 
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